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FOREWORD 

This  report  covers  research  performed  under  Air  Force  contract  F04611-82- 
C-0004  during  the  period  November  15,  1981  to  February  15,  1984.  The  high 
pressure,  high  temperature  CARS  spectrum  of  carbon  monoxide  was  investigated 
for  the  end  purpose  of  permitting  accurate  temperature  and  concentration 
measurements  of  this  species  in  combustion  media.  This  was  accomplished  by 
the  acquisition  of  high  quality  CARS  spectra  under  controlled  laboratory 
conditions  of  temperature  and  pressure  which  were  used  to  refine  and  verify  a 
theoretical  model  that  predicts  the  spectra  through  computer  simulation.  The 
computer  model  can  then  be  used  in  regression  analyses  on  experimental  spectra 
to  measure  accurately  the  temperature  and  concentration  of  the  combustion 
species. 

The  personnel  involved  with  this  contract  have  had  responsibilities  for 
different  phases  of  work.  John  H.  Stufflebeam  performed  the  experimental 
investigations  and  acquired  the  CO  CARS  spectra;  Robert  J.  Hall  was 
responsible  for  theoretical  aspects  of  high  pressure  CO  CARS  spectroscopy, 
including  the  computer  analyses  employed  to  model  the  expermintal  data.  Alan 
C.  Eckbreth  provided  the  overall  direction  of  the  experimental  and  theoretical 
aspects  of  this  work  and,  in  addition,  supplied  specific  solutions  to 
experimental  problems  that  were  encountered.  This  contract  was  monitored  by 
David  P.  Weaver  of  AFRFL. 
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SUMMARY 


Under  Contract  F0A61 1-82-C-0004  sponsored  by  the  Air  Force  Rocket 
Propulsion  Laboratory  (RPL),  the  United  Technologies  Research  Center  (UTRC) 
has  conducte^t^pectroscopic  investigation^  into  the  coherent  anti-Stokes 
Raman  spectroscopy  (CARS)  spectrum  of  CO.  CARS  is  a  remote  laser  diagnos¬ 
tic  technique  for  temperature  and  species  measurements  in  hostile  combus¬ 
tion  environments.  As  such  it  possesses  considerable  relevance  to  the  Air 
Force  in  the  general  areas  of  ballistics  and  propulsion.  This  final  report 
describes  the  results  of  investigations  conducted  into  the  effects  of  high 
pressure,  specifically  the  phenomenon  of  collisional  narrowing,  on  CARS 
spectra  from  which  temperature  and  density  information  derive.  Experimen¬ 
tal  studies  of  CARS  spectra  were  conducted  in  carbon  monoxide  in  a  heated, 
high  pressure  cell  over  the  range  300K-1500K,  1-100  atm.  The  spectra  were 
used  to  evaluate  the  UTRC  CARS  computer  code  supplied  to  RPL  for  the  calcu¬ 
lation  of  CO  CARS  frequency  distributions. 
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INTRODUCTION 

High  pressure  combustion  is  extremely  Important  in  a  variety  of  prac¬ 
tical  applications  of  Air  Force  relevance  such  as  propulsion,  airbreathing 
and  otherwise,  and  ballistics.  Experimental  diagnostics  of  high  pressure 
combustion  phenomena  are  Important  to  gain  the  understanding  required  to 
Improve  and  control  these  processes.  Very  little  is  known  about  high  pres¬ 
sure  combustion  processes  due  to  the  lack  of  suitable  diagnostic  tech¬ 
niques.  A  vast  archive  of  data  exists  from  insertion  of  thermocouples, 
pressure  transducers,  and  gas  sampling  probes  into  combustion  media.  How¬ 
ever,  physical  probes  are  not  only  confronted  with  survival  considerations, 
they  are  of  questionable  utility  since  their  presence  may  seriously  perturb 
the  phenomena  under  study.  Optical  techniques  appear  Ideally  suited  to 
diagnosing  such  phenomena,  but  many  are  inappropriate  due  to  poor  spatial 
resolution,  weak  signal  strength,  etc.  One  optical  technique,  coherent 
anti-Stokes  Raman  spectroscopy  or  CARS,  shows  considerable  promise  for  di¬ 
agnosing  high  pressure  combustion  phenomena  (Refs.  1,2).  CARS  is  spatial¬ 
ly  precise  (Ref.  3)  and  its  coherent  or  beamlike  nature  is  quite  amenable 
to  the  limited  optical  apertures  typical  of  high  pressure  combustion  cham¬ 
bers.  Furthermore,  CARS  exhibits  a  nonlinear  dependence  on  molecular  num¬ 
ber  density  and  increases  rapidly  in  signal  strength  as  the  gas  density  is 
elevated.  CARS  has  already  been  shown  to  be  practically  applicable,  having 
been  demonstrated  in  gas  turbine  combustors,  internal  combustion  engines, 
propellant  flames,  furnaces,  and  shock  tubes  (Ref.  4  and  references  con¬ 
tained  therein). 

In  CARS,  incident  laser  beams  at  frequencies  (Oj  and  (1)2  (termed  the 
pump  and  Stokes,  respectively),  interact  through  the  third-order,  nonlinear 
electric  susceptibility,  X'^^,  to  produce  the  coherent  CARS  radiation  at  fre¬ 
quency  Uj  “  2a)j  -  (1)2.  If  the  frequency  difference  u)j  -  u)2  coincides 
with  a  Raman-active  vibrational  or  rotational  mode  of  a  certain  species, 
then  the  CARS  radiation  is  resonantly  enhanced  and  uniquely  characteristic 
of  that  species.  In  CARS,  measurements  of  temperature  are  performed  from 
the  spectral  shape,  i.e.,  the  CARS  intensity  distribution  with  frequency 
(Ref.  5).  Species  density  can  be  determined  from  the  spectral  shape  as 
well  in  certain  concentration  ranges,  but,  in  general,  density  measurements 
are  based  upon  the  strength  of  the  CARS  signal  (Ref.  6).  The  intensity  of 
CARS,  however,  is  quite  linewldth  dependent.  Collision  processes  alter  the 
spectra  by  perturbing  the  rovibronic  linewidths.  As  the  gas  density  in¬ 
creases,  a  variety  of  linewldth  phenomena  occur  such  as  pressure  broaden¬ 
ing,  shifting,  motional  or  Dicke  narrowing,  and  collisional  narrowing  which 
complicate  the  calculation  of  the  CARS  spectrum.  For  CARS,  the  collisional 
narrowing  phenomenon  is  important  for  pressures  above  one  atmosphere.  Col¬ 
lisional  narrowing  occurs  when  adjacent  Raman  transitions  have  been  pres¬ 
sure  broadened  to  the  extent  that  they  overlap.  The  overlap  allows  commu¬ 
nication  between  lines  and  they  subsequently  coalesce  or  collapse  to  a  nar¬ 
rower  bandwidth.  Accurate  modeling  of  high-pressure,  high  temperature  CARS 
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spectra  Is  necessary  to  extract  temperature  and  density  information  and, 
thus,  is  fundamental  to  the  deployment  of  the  technique  as  a  diagnostic. 
Well  controlled,  laboratory  experiments  provide  numerical  input  to  the  mod¬ 
el. 


Under  Air  Force  sponsorship.  Contract  F0461 1-82-C-0004,  UTRC  has  in¬ 
vestigated  the  phenomenon  of  pressure-induced,  l.e.  colllsional,  narrowing 
in  carbon  monoxide.  CO  CARS  spectra  have  been  obtained  over  the  range 
1-100  atm.,  300K-1500K.  Specific  combinations  of  these  parameters  were 
selected  for  the  acquisition  of  CARS  spectra  as  Indicated  in  Fig.  1.  Ex¬ 
cellent  agreement  in  CO  has  been  obtained  between  the  experimental  spectra 
and  those  predicted  by  the  colllsional  narrowing  theory. 

The  next  section  is  a  summary  of  the  theoretical  model  of  colllsional 
narrowing  and  its  application  to  the  prediction  of  CO  CARS  spectra.  Sub¬ 
sequent  to  that,  the  experimental  facility  is  described,  together  with  the 
techniques  employed  to  obtain  the  high  quality  CO  spectra.  The  results  of 
this  investigation  are  then  presented  and  comparisons  are  made  between  the 
predictions  from  the  theoretical  model  and  the  experimental  spectra.  Dis¬ 
cussions  of  these  results  are  Included  with  specific  reference  to  valida¬ 
tion  of  the  computer  code  for  CO.  Suggestions  for  future  work  are  outlined 
and  finally,,  the  experimental  data  obtained  during  this  Investigation  are 
contained  as  an  attached  appendix. 
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COLLISIONAL  HARROWING  OF  CARS  SPECTRA 

The  basic  theory  of  CARS  for  the  case  in  which  none  of  the  field  fre¬ 
quencies  encroaches  upon  electronic  resonances  in  the  probed  molecule  has 
been  given  in  many  publications  (Refs.  1,2, 7-9).  It  is  well  known  that  the 
anti-Stokes  frequency  is  proportional  to  the  squared  modulus  of  the  third 
order  electric  susceptibility.  Since  CARS  temperature  measurements  and,  in 
certain  ranges,  concentration  measurements  are  performed  from  the  shape  of 
the  signature  (i.e.  the  variation  of  intensity  with  frequency),  the  propor¬ 
tionality  constants  are  not  of  general  interest.  The  third-order  suscepti¬ 
bility  can  in  turn  be  represented  as  a  sum  of  frequency-dependent,  complex 
terms  associated  with  Raman  resonances  plus  a  nonresonant  background  con¬ 
tribution  that  will  be  essentially  real  and  dispersionless  if  one  and  two 
photon  resonances  are  avoided  in  all  species  in  the  mixture  of  Interest. 

Of  particular  theoretical  interest  for  this  contract  is  the  effect  of  line 
overlap  at  high  pressures,  which  gives  rise  to  the  phenomenon  of  collislon- 
al  narrowing.  This  is  distinct  from  Dlcke,  or  motional,  narrowing  wherein 
discrete  thermal  velocity  groups  blend  together  through  elastic,  momentum- 
transferring  collisions.  A  CO  vibrational  Q-branch  consists  of  many  close¬ 
ly  spaced  transitions,  corresponding  to  each  rotational  quantum  number  j 
that  are  split  very  slightly  by  vibration-rotation  Interaction.  Different 
vibrational  bands  in  CO  will  have  a  large  splitting  due  to  vibrational  an- 
harmonicity,  but  within  each  band  the  overlap  can  become  very  strong  at 
modest  pressures.  Under  such  conditions,  it  no  longer  suffices  to  repre¬ 
sent  the  resonant  part  of  the  susceptibility  as  a  simple  sum  of  contribu¬ 
tions  from  Independently  broadened  transitions,  and  very  significant  errors 
will  be  incurred  if  the  effects  associated  with  colllslonal  narrowing  are 
Ignored . 


CARS  Colllslonal  Narrowing  Formalism 

Generalized  treatments  of  the  colllslonal  narrowing  effect  have  been 
given  by  several  authors  (Refs.  10-12),  but  the  analyses  that  are  of  most 
relevance  to  the  present  work  are  the  extension  to  spontaneous  Raman  Q- 
branches  by  Alekseyev  and  Sobelmann  (Ref.  13),  and  the  extension  to  CARS 
spectroscopy  of  vibrational  Q-branches  by  Hall,  et.al.  (Ref.  lA).  In  the 
latter  formulation,  the  resonant  Raman  contibution  to  the  third  order  sus¬ 
ceptibility  is  given  in  the  following  form: 


(1) 


where  the  k  and  j  summations  are  over  all  Raman  transitions;  N  Is  the 
probed  s|)ecies  number  density ,  a  denotes  polarizability  matrix  element, 
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(o ) 

Ap.j  denotes  the  unperturbed  population  difference  between  the  initial 
and  final  states  in  each  transition,  and  the  elements  of  the  matrix  C, 
which  has  to  be  inverted,  are  given  by: 

(-^  -'^k)  +  Xkj  (2) 

where  u^,  co2>  and  are  the  frequencies  of  the  pump  laser,  the  Stokes 
laser,  and  Raman  transition  Q(k),  respectively;  is  the  Raman  linewidth 
of  transition  Q(k);  Aj^  is  the  pressure-induced  shift  of  the  transition,  and 
denotes  an  off-diagonal  linewidth  element  which  plays  a  key  role  in  de¬ 
termining  the  extent  of  collislonal  narrowing. 

Normally  enough  Information  is  available  to  fill  out  the  diagonal  ele¬ 
ments  of  the  G-matrlx,  but  the  off-diagonal  elements  are  a  more  complex  mat¬ 
ter.  They  can  be  represented  in  terms  of  the  scattering  matrix  (S)  for  in¬ 
elastic  energy  transfer  (Refs.  13,14)  and  certain  approximations  can  be  in¬ 
troduced.  For  example,  if  vibrational  energy  transfer  rates  are  slov 
they  probably  are  in  CO,  it  can  be  shown  that  the  off-diagonal  eleme  only 

couple  transitions  with  the  same  initial  and  final  vibrational  state  Thus, 

the  G-matrix  will  be  block  diagonal,  with  one  block  for  each  vibrat:  1  band 

of  interest.  If  it  is  also  assumed  that  the  scattering  matrix  is  on  /  akly 
dependent  on  vibrational  state,  which  means  that  the  Raman  linewidths  are  not 
dependent  on  vibrational  band,  then  It  can  be  shown  that  the  off-diagonal  ele¬ 
ment  is  minus  the  first  order  rate  coefficient  for  Inelastic  rotational 
energy  transfer  in  the  direction  J  k  (Refs.  13,14).  Vibrational  dephaslng 
effects  can  arise  from  purely  vibrational  collislonal  processes  and  from  ro¬ 
tational  dephasing  collisions  if  the  effect  of  vibration-rotation  Interaction 
is  significant.  These  effects  will  give  rise  to  ordinary  pressure  broadening 
behavior  even  when  adjacent  lines  are  strongly  overlapped,  and  will  counteract 
the  narrowing  effect  of  the  rotationally  inelastic  collisions.  The  rotation- 
ally  Inelastic  collisions  always  contribute  to  the  Isolated  linewidth,  but 
when  the  pressure  has  In^'raased  to  the  point  where  adjacent  lines  are  over¬ 
lapped,  they  give  rise  to  a  narrowing  of  the  spectrum.  All  the  evidence  is 
that  the  pure  vibrational  processes  make  small  contributions  to  the  isolated 
linewidths,  however. 

Thus,  the  problem  of  specifying  the  linewidth  parameters  which  govern  the 
narrowing  process  reduces,  to  a  good  approximation,  to  one  of  specifying  the 
rotational  relaxation  rate  matrix;  this  means  that  certain  relationships  such 
as  detailed  balance  and  S-matrix  unitarity  can  be  invoked  to  reduce  the  number 
of  unknowns : 
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AO)  (O) 

^hk^jk  ”  ^jj  ^kj 


(4) 


However,  these  relationships  still  do  not  specify  the  relaxation  rate  matrix 
entirely,  even  If  one  knows  the  diagonal,  isolated  llnewldths.  In  principle, 
ab  Initio  scattering  calculations  could  be  employed  to  generate  this  matrix, 
but  there  are  many  problems  with  such  an  approach.  Accurate  trajectory  cal¬ 
culations  require  good  Intermolecular  potential  surfaces,  and  these  are  well 
known  only  for  the  most  simple  systems.  Also,  such  calculations  would  be 
very  time-consuming,  and  would  have  to  be  performed  at  each  temperature  and 
composition  of  Interest.  Even  in  the  event  that  this  could  be  done.  It  Is 
unlikely  that  the  results  could  be  expressed  In  a  succinct  functional  form 
that  would  facilitate  their  application  to  practical  systems.  For  CO  these 
rates  could  In  principle  be  measured  experimentally,  but  again,  only  a  few 
such  state-to-state  measurements  are  available.  The  approach  adopted  In  this 
program  has  been  to  exploit  the  considerable  success  enjoyed  by  simple,  seml- 
emplrlcal  scaling  laws  In  correlating  rotational  relaxation  rate  data.  A  good 
review  of  this  field  has  been  given  by  Brunner  and  Pritchard  (Ref.  15).  Typi¬ 
cally,  the  state-to-state  rates  are  given  by  a  simple  scaling  law  which  is  a 
function  of  the  translational  energy  defect  associated  with  the  relaxation 
process.  The  most  well-known  laws  are  the  exponential  gap  law  (Ref.  16)  and 
the  inverse  power  law  (Ref.  17).  The  scaling  laws  investigated  for  CO  will  be 
described  In  detail  in  the  following  section;  next  however,  some  simplified 
calculatlonal  algorithms  will  be  discussed. 

Confirmation  of  the  basic  correctness  of  the  G-matrlx  approach  has  been 
given  mainly  by  comparisons  with  experimental  N2  CARS  spectra  for  pressures 
ranging  up  to  100  atmospheres  (Refs.  14,17).  The  work  of  Rosasco  and  co¬ 
workers  (Ref.  18),  who  used  the  same  approach  in  low  pressure  N2  and  CO  in¬ 
verse  Raman  spectroscopy  (IRS)  also  serves  to  confirm  the  correctness  of  the 
basic  theory.  One  calculatlonal  difficulty  associated  with  modelling  high 
pressure  CARS  spectra  is  that  the  size  of  the  G-matrlx  can  become  quite  large, 
particularly  at  flame  temperatures;  because  this  matrix  must  be  Inverted  at 
each  frequency  mesh  point  one  selects  to  perform  a  calculation,  the  computer 
run  times  can  become  extremely  long.  It  is  desirable,  therefore,  to  have  sim¬ 
pler  calculatlonal  algorithms.  One  approach  useful  at  moderate  pressures  when 
the  Isolated  line  approximation  is  no  longer  valid  is  one  based  upon  a  pertur¬ 
bation  expansion  of  the  G-matrix  as  a  function  of  pressure.  In  Ref.  14  it  was 
shown  how  the  Inverse  matrix  could  be  so  expanded,  and  the  result  for  the 
first  order  correction  to  the  pressure  is: 
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This  approach  was  also  used  in  Ref.  (18).  If  one  retains  the  second  order 
term,  then  it  is  possible  to  employ  the  following  second  order  perturbation 
expansion: 
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It  has  been  assumed  in  Eqs.  8,9  that  (Oj  -  u)2  “  for  any  corrective 
terras  to  to  be  important.  This  perturbation  expansion  is  particularly  use¬ 
ful  for  moderate  pressures  in  hot  gases,  but  is  not  applicable  over  the 
whole  range  of  conditions  of  interest  in  this  study.  It  has  the  added  ad¬ 
vantage  that  it  is  calculationally  very  efficient. 


Another  efficient  colllslonal  narrowing  algorithm  that  has  been  found 
to  give  reasonable  results  in  N2  is  associated  with  Gordon's  extended  J- 
dif fusion  theory  of  rotational  motion  (Refs.  19-21).  This  simplified  theo¬ 
ry  assumes  strong  collisions  in  which  the  magnitude  and  orientation  of  the 
angular  momentum  are  randomized;  there  are  thus  no  selection  rules  for  ro¬ 
tational  relaxation,  and  the  probability  of  transfer  from  rotational  state 
1  to  rotational  state  J  is  just  proportional  to  the  Boltzmann  population  of 
state  j.  In  formulating  the  third  order  susceptibility  for  this  model,  it 
is  possible  to  circumvent  the  matrix  inversion  step  (Refs.  20-21),  and  the 
result  is  the  very  simple  expression  given  by  the  following: 
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The  effect  of  the  narrowing  is  controlled  by  the  simple  summation  term  in 
the  denominator.  In  addition  to  having  some  success  with  N2,  the  rotational 
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diffusion  theory  has  been  successfully  applied  by  Papineau  and  Pealat  (Ref.  22) 
to  modelling  CO2  vibrational  CARS  spectra  at  room  temperature;  because  the 
vibration-rotation  Interaction  is  so  small  in  this  molecule,  substantial 
collislonal  narrowing  effects  are  encountered  even  at  one  atmosphere.  The 
great  advantage  of  this  expression  relative  to  the  full  G-matrix  algorithm 
is  its  efficiency;  it  requires  essentially  no  more  run  time  than  an  iso¬ 
lated  line  calculation. 

The  UTRC  high  pressure  code  has  as  options  these  various  approximations 
to  the  line  overlap  problem.  .t,e  user  may  specify  that  calculations  be 
based  upon  either  the  Isolated  line  approximation,  the  small  perturbation 
expression,  the  Gordon  rotational  diffusion  calculation,  or  the  full  G- 
matrlx  inversion. 

CO  Llnewidth  Parameters 

The  most  fundamental  parameters  needed  for  modelling  CARS  spectra  are 
the  Isolated  linewidths,  Fj.  For  CO,  the  Q-branch  Raman  linewidths  are 
now  known  with  great  precision  from  the  Inverse  Raman  spectroscopy  measure¬ 
ments  of  Rosasco  and  co-workers  (Ref. 18)  and  are  shown  in  Fig.  2.  These 
linewidths  are  in  very  good  agreement  with  infrared  P-  and  R-branch  mea¬ 
surements  by  BelBruno,  et.  al.  (Ref.  23);  because  the  IRS  measurements  did 
not  extend  beyond  j=20,  the  BelBruno  P-branch  linewidths  were  employed  out 
to  j“25  to  give  a  set  that  could  be  used  for  300K  CARS  calculations.  BelBruno, 
et.  al.  also  measured  the  CO  vibrational  overtone  widths,  and  found  only  a 
slight  temperature  dependence;  all  of  these  results  are  an  indication  that  it 
is  mainly  rotational  relaxation  that  makes  the  dominant  contribution  to  the 
linewidths.  The  slight  band  dependence  of  the  overtone  measurements  may 
suggest  a  small  elastic  or  vibrational  contribution.  The  result  that  the 
Raman  and  infrared  linewidths  agree  well  is  also  potentially  useful,  be¬ 
cause  it  indicates  that  Infrared  llnewidth  data  concerning  temperature  and 
concentration  dependences  could  be  useful  to  fill  in  gaps  in  the  Raman 
data. 


Unfortunately,  there  is  very  little  data  on  the  temperature  dependence  i 
of  any  of  the  pure  CO  linewidths.  In  general,  one  can  express  the  tempera-  ) 
ture  dependence  in  the  form:  I 


(11) 

where  the  temperature  exponent  hjis  expected  to  be  j-dependent.  Sun  and 
Griffiths  (Ref.2A)  state  that  rij  may  take  on  values  in  the  range  0.5  to 
1.0,  which  presumably  reflects  the  range  expected  as  the  interaction  varies 
from  long-range  multipole  to  short-range,  hard-sphere.  Sun  and  Griffiths 
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report  that  their  measurements  suggest  an  exponent  of  approximately  0.66 
for  lines  P(25)  and  P(31)  over  the  temperature  range  from  ambient  to  180  C, 
as  opposed  to  the  previously  accepted  value  of  approximately  0.75.  How¬ 
ever,  this  represents  a  relatively  small  temperature  range  for  rotational 
levels  that  are  not  of  great  significance  for  low  temperature  CARS  calcula¬ 
tions.  Varghese  and  Hanson  (Refs.  25,26)  have  investigated  the  effect  of 
collision  partner  on  room  temperature  CO  P-branch  llnewidths,  and  from  mea¬ 
surements  of  these  widths  in  a  flame  have  made  some  deductions  about  the 
temperature  exponent.  For  CO  broadened  by  N2,  the  exponent  ranges  from 
roughly  0.8  at  low  j  to  0.5  at  high  j.  However,  none  of  this  gives  one  in¬ 
formation  directly  about  the  temperature  dependences  of  the  pure  CO  P- 
branch  widths.  It  also  seems  difficult  to  resort  to  theoretical  arguments 
to  deduce  how  the  pure  broadening  exponents  might  differ  from  those  for  N2. 

Ab  initio  llnewidth  calculations  for  pure  CO  similar  to  those  performed 
successfully  for  pure  N2  using  the  theory  of  Bonamy  and  Robert  (Ref.  27) 
would  be  a  very  useful  guide,  but  are  very  time-consuming  and  were  judged 
to  be  outside  the  scope  of  this  contract.  A  compromise  was  made  in  which 
use  was  made  of  the  theoretical  calculations  for  pure  N2  Q-branch  line- 
widths,  in  which  to  a  good  approximation 


The  parameters  A  and  B  were  determined  from  a  least  squares  fit  to  the  room 
temperature  data  in  Fig.  2,  and  the  temperature  exponents  and 
were  assumed,  lacking  better  information,  to  be  equal  to  those  for  pure  N2 
(0.7  and  1.45).  The  exponent  is  certainly  not  Inconsistent  with  the 
information  previously  presented.  The  low  j  lines  will  have  the  strongest 
temperature  dependence,  and  the  high  j  lines  will  vary  more  slowly  due  to 
the  fact  that  the  slope  decreases  rapidly  with  temperature.  Generally,  the 
greatest  sensitivity  to  llnewidth  occurs  at  low  pressures. 

With  reasonable  estimates  for  the  diagonal  llnewidths,  the  off-diag¬ 
onal  elements  can  then  be  calculated  if  one  assumes  a  particular  functional 
form  for  the  state-to-state  rate  matrix.  It  is  likely  that  each  molecule 
will  represent  a  different  situation,  with  no  one  rate  model  having  univer¬ 
sal  applicability.  In  a  semi-empirical  fashion,  the  particular  scaling  or 
rate  law  which  best  satisfies  all  of  the  data  on  the  pressure  dependence 
must  be  found.  There  is  also  a  choice  of  ways  to  fill  out  the  rate  matrix 
for  a  given  functional  form.  One  alternative  is  to  assume  the  functional 
form  for  endothermic  rates  (j  -  k,  k>j),  and  then  use  detailed  balance  and 
the  conservation  of  probability  relationships  to  deduce  exothermic  rates 
(j  ^  k,  k<j),  in  this  way  working  down  the  rows  of  the  rate  matrix.  The 
method  pursued  here  has  been  to  least  squares  fit  a  rate  coefficient  model 
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satisfying  detailed  balance  in  advance  to  the  knovm  isolated  linewidths 
through  the  relationship: 


where  the  theoretical  linewidths  are  given  by 


1 

2  i 


Xij 


(13) 


(14) 


The  Isolated  linewidths  used  in  the  calculation  must  then  be  identified 
with  the  theoretically-deduced  values;  usually  there  will  be  no  problem 
with  this  because  the  fits  are  quite  good. 

For  N2,  all  of  the  pressure  dependence  data  presently  available  In¬ 
dicates  that  a  simple,  population-weighted  Inverse  power  law  of  the  form 

=  (15) 

does  a  good  job.  The  fits  to  the  isolated  linewidths  are  good,  and  the 
latter  can  be  regarded  as  reasonably  well  known  for  N2.  This  power  law, 
hereafter  referred  to  by  the  acronym  SPL,  gives  good  agreement  with  the 
IRS  measurements  of  Ref.  18  at  room  temperature  and  modest  pressures  and 
also  with  UTRC  spectra  taken  at  100  atmospheres  and  a  temperature  of  ap¬ 
proximately  1700K  (Ref.  17).  This  model  was  also  investigated  for  CO 
and,  as  will  be  seen,  gives  generally  good  agreement  with  the  experiments. 
Another  rate  model  that  was  originally  proposed  by  Dexhelmer  (Ref.  28)  was 
also  investigated  as  part  of  this  contract: 
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where  the  factor  b  is  now  an  additional  fit  parameter.  This  is  termed  an 
exponential  po%rer  law  model,  and  will  hereafter  be  referred  to  as  the  EPL. 

The  fits  to  the  isolated  liomwldths  obtained  from  these  two  models  are 
also  exhibited  in  Fig.  2.  The  SPL  gives  a  generally  good  fit,  with  devia¬ 
tions  from  experiment  most  pronounced  at  low  J;  the  disagreeoient  here  is 
not  necessarily  the  result  of  a  defect  in  the  model  because  elastic  contri 
buttons  could  be  important  at  the  smaller  J  values.  The  fit  is  similar  to 
that  obtained  by  BelBruno,  et.al.  (Ref.  23)  who  inverted  their  P-branch 
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llnewldths  to  obtain  rate  coefficients  from  a  much  more  elaborate  rate  mod¬ 
el.  With  the  extra  fitting  parameter,  the  EPL  gives  a  fit  that  is  essen¬ 
tially  perfect,  as  shown  in  Fig.  2.  This  does  not  mean  that  it  is  a  bette 
representation  of  the  relaxation  matrix,  however.  As  will  be  seen,  it  is 
difficult  to  choose  between  the  two  models  based  on  the  high  pressure  data 
It  should  be  noted  that  the  EPL  does  not  appear  to  be  consistent  with  the 
UTRC  high  temperature,  high  pressure  data  of  Ref.  17. 
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EXPERIMENTAL  APPARATUS  AND  TECHNIQUES 
CARS  Apparatus 

A  schematic  of  the  apparatus  used  to  obtain  the  CARS  spectra  Is  shown 
In  Fig.  3.  The  output  of  a  Quanta  Ray  Nd:YAG  laser  Is  frequency  doubled  to 
generate  a  horizontally  polarized,  "primary”  pump  beam  at  5320A  The 

lOHz  laser  has  an  average  output  of  2.2  Watts  at  5320A  with  a  pulse  dura¬ 
tion  of  10  nanoseconds.  A  beamsplitter  separates  30Z  from  u ^  to  pump  the 
broadband  Stokes  dye  laser  oscillator  and  amplifier.  Galilean  telescopes 
are  provided  to  control  beam  waists  and  the  focal  zone  locations  of  and 
(0  2  which  are  combined  on  the  dlchroic  mirror  and  focused  Inside  the  high 
pressure  test  cell.  CARS,  at  (03,  is  generated  In  the  focal  volume,  and 
all  three  frequencies  (u  ^  ,  (1)2,  )  are  recollimated  after  exiting  the 

cell.  A  second  dlchroic  mirror  separates  (1)3  from  ti>  2  before  the  sig¬ 

nal  Is  incident  on  the  slit  of  a  double,  one-meter  monochromator  with  a 
limiting  resolution  of  ~0.3  cm  .  The  m  and  u  2  beams  are  trapped  after 
reflection  from  the  second  dlchroic.  A  He-Ne  laser  is  shown  whose  output 
Is  coincident  with  the  (Dj^,  u  2  path  and  used  for  optical  alignment.  A 
beamsplitter  In  the  leg  provides  a  reference  signal  to  a  photomulti¬ 
plier  which  Is  used  to  normalize  the  spectrally  dispersed  CARS  signal  from 
the  monochromator.  The  two  signals  are  ratioed  in  a  boxcar  averager  to 
account  for  variations  In  laser  Intensity  during  the  scan  of  the  CARS  spec¬ 
trum.  At  high  gas  pressures  It  Is  necessary  to  attenuate  the  lasers  to 
avoid  gas  breakdown  and  stimulated  Raman  gain  which  would  Interfere  with 
and  perturb  the  CARS  signal.  The  attenuation  is  achieved  through  a  Fresnel 
reflection  (~  1.4%)  from  a  flat  surface.  A  prism  is  used  for  this  purpose, 
as  shown  In  Fig.  3,  to  avoid  the  back  surface  reflection  that  would  result 
from  a  plane  parallel  element.  Alternate  methods  of  attenuation  such  as 
absorbing  neutral  density  filters  or  reduction  In  the  flashlamp  energy  in¬ 
put  to  the  NdrYAG  laser  are  less  desirable.  They  cause  changes  In  the 
focal  zone  location  due  to  thermal  leasing  In  the  filters  or  a  change  In 
divergence  of  the  output  of  the  Nd:YAG  laser,  and  the  focal  shift  can  be 
great  enough  to  cause  damage  to  the  sapphire  rods.  A  Fabry-Perot  etalon  Is 
Included  to  monitor  the  llnewldth  of  the  primary  pump  beam.  The  convolu¬ 
tion  of  and  U2  determine  the  resolution  of  the  CARS  spectra  at  U3. 

An  angle  tuned  etalon  option  on  the  Nd:YAG  laser  results  in  an  line- 
width  of  ~0.47  cm”^  and  the  external  etalon  is  used  to  detect  any  shift 
from  this  value  which  would  degrade  the  resolution  of  the  CARS  spectra.  A 
vldicon  detector  is  mounted  on  the  spectrograph,  used  in  the  single  mono¬ 
chromator  configuration,  with  an  overall  instrument  resolution  of  1.1  cm 
which  Is  adequate  for  broadband  spectra.  This  detector  can  cover  a  broad 
spectral  range  (>200  cm~^ )  and  Is  used  to  monitor  the  broadband  Stokes 
dye  laser  spectral  position  for  the  CO  experiments.  The  Instrument  is  In¬ 
terfaced  to  a  PDF  11/35  computer  to  control  Its  operation  and  store  the 
spectral  data  for  subsequent  digital  analysis.  A  photograph  of  the  UTRC 
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FIG.  4 


HIGH  PRESSURE  CARS  FACILITY 
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High  Pressure  Facility  Is  shown  In  Fig.  A  with  many  of  the  components 
Identified  that  were  described  in  the  previous  discussion. 

High  Pressure  Cell 

To  provide  well  controlled  conditions  of  temp>erature  and  pressure  for 
the  acquisition  of  CO  CARS  spectra,  the  high  pressure  vessel  schematically 
depicted  In  Fig.  5  was  utilized.  It  is  a  static.  Internally  heated  cell 
that  is  rated  for  temperatures  to  1750K  and  pressures  to  5000  pslg.  It  is 
equipped  with  a  pressure  gauge,  type  S  thermocouples,  and  an  optical  pyro¬ 
meter  to  measure  the  pressure  and  temperature  of  the  sample.  Only  the  cen¬ 
tral  15  cm  is  heated  and  sapphire  rods  extend  into  the  heated  region  to  re¬ 
duce  thermally-induced  density  gradient  effects  in  the  gas.  The  stepped, 
hollow  ceramic  spacer.  Identified  in  the  figure,  is  used  as  a  target  for  an 
optical  pyrometer  to  measure  the  temperature  of  the  wall  adjacent  to  the 
gas.  The  pyrometer  reading  is  analytically  corrected  for  the  emisslvity  of 
the  alumina  and  transmission  through  the  sapphire  and  quartz.  In  the  re¬ 
gion  where  measurements  by  the  thermocouple  and  pyrometer  overlap,  their 
temperature  data  exhibit  a  linear  relationship.  This  linear  relationship 
is  extended  to  provide  temperature  data  from  the  thermocouple  in  the  region 
where  pyrometer  measurements  are  Impossible,  i.e.  T  <  950K. 

The  axial  temperature  profile  of  the  cell  has  been  measured  at  atmo¬ 
spheric  pressure  and  the  results  are  shown  in  Fig.  6.  The  profiles  for  two 
different  central  temperatures  are  shown.  As  expected,  the  temperature 
peaks  In  the  center  of  the  cell  and  smoothly  decays  to  room  temperature  at 
the  Input  window.  Also  indicated  in  Fig.  6  is  the  CARS  interaction  length, 
70  mm,  measured  for  the  two  beam  BOXCARS  geometry  described  in  the  next 
paragraph.  The  data  show  that  the  temperature  varies  only  ±  3.6X  over  the 
CARS  interaction  length. 

A  peculiarity  of  the  carbon  monoxide  molecule  was  encountered  during 
this  Investigation.  It  Involved  the  disproportionation  of  CO: 

2C0  -  CO2  +  C  (17) 

Ref.  29  indicates  that  this  reaction  is  thermodynamically  favored  for  tem¬ 
peratures  between  700K  and  lOOOK,  in  the  presence  of  an  active  surface. 

The  disproportionation  was  observed  during  the  CO  investigations,  over  the 
temperature  range  Indicated  and  only  when  the  sapphire  rode  were  in  the 
cell.  This  initially  precluded  the  acquisition  of  CO  CARS  spectra  above 
900K  because  the  CO  decomposed  to  CO2  and  carbon.  The  carbon  coated  out 
on  the  faces  of  the  rods  which  caused  absorption  and  attenuation  of  the 
pump  lasers  so  that  no  CARS  could  be  generated.  This  result  was  verified 
by  SIMS  (Secondary  Ion  Mass  Spectroscopy)  analysis  of  the  deposits  on  the 
sapphire  rods.  After  considerable  experimental  effort  and  collaboration, 
the  simple  solution  was  to  heat  the  cell  to  the  desired  temperature  at  1  atm. 


INTERNALLY  HEATED  PRESSURE  VESSEL  FOR  CARS  DIAGNOSTICS 
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pressure,  then  fill  the  hot  cell  to  the  pressure  desired  for  the  CARS 
spectra  and  wait  —30  minutes  for  the  gas  to  equilibrate.  Some  dispropor¬ 
tionation  still  occurred  but  not  at  the  spatial  location  of  the  CARS  sample 
volume.  The  center  of  the  cell  was  above  the  disproportionation  tempera¬ 
ture  so  the  ends  of  the  rods  remained  relatively  clear  which  allowed  ac¬ 
quisition  of  the  high  pressure,  high  temperature  spectra.  Inspection  of 
the  rods  after  this  exercise  revealed  a  band  of  soot  approximately  3  cm 
from  the  end  of  the  sapphire  rods  -  apparently  the  location  of  the  optimum 
temperature  for  disproportionation. 

Phase  Matching  Approaches 

The  limited  optical  access  provided  by  the  test  cell  dictates  either  a 
coaxial  alignment,  or  tightly  packed  BOXCA^-S  geometry,  of  the  pump  and 
Stokes  lasers  used  to  generate  CARS.  Acquisition  of  CO  spectra  at  room 
temperature  or  at  elevated  temperatures  above  10  atm.  pressure  employed  a 
two  beam,  three  dimensional  phase-matching  configuration  called  USED  CARS 
(Ref.  30).  In  this  two  beam  approach,  Che  m  2  Stokes  beam  is  coaxial  with 
and  inside  the  annular  u  ^  pump  beam.  The  annular  pump  results  from  Che 
unstable  resonator  employed  on  Che  Nd:YAG  laser.  The  pump  and  Stokes  beams 
do  not  overlap  until  the  focal  region  where  the  pump  beam  transforms  or 
fills  in.  The  resulting  resolution  elememnt  is  a  cylinder  of  0.1  mm  diame¬ 
ter  and  100  mm  long.  For  295K  spectra  there  are  no  thermal  gradients  to 
complicate  interpretation  of  the  spectra  and  this  resolution  is  adequate. 

In  the  high  pressure,  high  temperature  experiments,  (P  >  10  atm.),  thermal 
gradients  cause  refractive  index  changes  in  the  gas  which  steer  the  laser 
beams  to  the  extent  they  cannot  be  propagated  through  the  cell.  Turbulence 
is  also  created  in  the  gas  and  destroys  the  spatial  coherence  of  the  pump 
beams  which  inhibits  the  generation  of  CARS.  These  effects  are  reduced 
significantly  by  the  insertion  of  sapphire  rods  and  high  quality  CARS 
spectra  are  observed.  The  presence  of  the  sapphire  rods  also  improves  the 
axial  resolution  by  limiting  Che  gas  sample  to  51  mm  in  length.  The  sap¬ 
phire  rods  were  not  used  for  the  spectra  at  pressures  from  1  to  10  atm. 
and  elevated  temperatures  because  pump  beam  intensities  that  would  damage 
the  rods  were  required  to  obtain  high  quality  spectra.  In  this  case  the 
best  axial  resolution  is  achieved  with  a  two  beam  BOXCARS  arrangement 
(Ref.  31)  in  which  the  Stokes  beam  is  located  on  the  circumference  at  one 
side  of  Che  annular  u ^  pump  beam.  The  measured  axial  interaction  length 
for  this  case  is  70  mm. 

The  two  beam  BOXCARS  geometry  is  Implemented  with  an  optical  flat, 
shown  in  Fig.  3,  Chat  is  inserted  after  the  Stokes  beam  amplifier  to  dis¬ 
place  (1)2  to  the  edge  of  (02  .  The  BOXCARS  formed  in  this  geometry  is  spa¬ 
tially  filtered  from  the  pump  and  Stokes  beams  as  well  as  from  any  collln- 
early  phase  matched  contributions  to  m j  since  the  colllnear  contributions 
will  have  poor  axial  resolution.  An  added  benefit  in  this  geometry  is  that 
less  (I)  2  is  incident  on  the  slit  of  the  monochromator  because  of  the  spa¬ 
tial  filtering.  This  reduces  the  noise  component  of  the  monochromator  out¬ 
put  and  results  in  higher  quality  spectra. 
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RESULTS  AND  DISCUSSION 
Experimental  Results 

CARS  spectra  of  carbon  monoxide  were  acquired  over  the  range  300K  - 
1500K,  1  -  100  atm.,  using  the  experimental  facility  and  techniques  de¬ 
scribed  above.  A  more  descriptive  synopsis  of  the  actual  parameters  inves¬ 
tigated  is  shown  in  Fig.  1  and  the  CO  spectra  themselves  are  Included  in 
the  appendix  of  this  report. 

Comparison  of  Theory  and  Experiment 

The  starting  point  for  modelling  the  experimental  spectra  is  the  295K 
spectrum  at  one  atmosphere  pressure,  as  given  in  Fig.  7.  It  is  apparent 
that  all  of  the  theoretical  approaches  reproduce  this  spectrum,  although 
the  EPL  does  have  a  tendency  to  predict  more  narrowing  than  is  observed. 

As  will  be  seen,  this  is  the  case  in  most  of  the  high  pressure  spectra. 

The  Gordon  model  gives  results  so  close  to  the  Isolated  line  and  SPL  pre¬ 
diction  that  it  is  not  shown  in  Fig.  7.  It  is  clear  that  there  is  no  mod¬ 
elling  problem  associated  with  the  295K  data  at  one  atmosphere. 

There  is  a  small  problem  associated  with  modelling  the  hot  CO  spectra 
at  one  atmosphere,  however,  with  diagreement  found  in  the  spectral  vicinity 
of  the  fundamental  bandhead.  An  example  is  shown  in  Fig.  8,  which  shows  the 
observed  pure  CO  spectrum  at  one  atmosphere  and  a  measured  temperature  of 
1265K.  The  somewhat  surprising  feature  of  the  spectrum  is  the  prominent 
feature  near  the  fundamental  bandhead;  the  rest  of  the  spectrum  shows  the 
expected  Q-branch  features  on  the  hot  side  of  the  fundamental,  as  well  as 
the  hot  band  with  no  unusual  features.  As  the  Appendix  shows,  this  band- 
head  shape  is  seen  in  all  the  hot  spectra.  Theoretical  calculations  using 
the  baseline  Raman  llnewidths  do  not  predict  that  this  feature  should  be  as 
significant  as  it  appears  to  be  experimentally,  although  there  is  a  small 
theoretical  peak  at  this  spectral  location.  Fig.  9.  Theoretically,  this 
peak  seems  to  arise  both  from  the  ordinary  cross-terms  interference  typical 
of  CARS  (that  is,  an  isolated  line  calculation  predicts  it),  and  from  col- 
lisional  narrowing  effects.  It  is  stronger,  when  as  in  Fig.  9,  collisional 
narrowing  is  Included  using  the  EPL  model.  If  the  Raman  llnewidths  in  the 
vicinity  of  the  bandhead  are  increased  by  assuming  a  temperature  exponent 
equal  to  0.5  instead  of  0.7,  the  theoretically  predicted  feature  does  be¬ 
come  more  prominent,  as  shown  in  Fig.  9,  but  is  still  in  disagreement  with 
experiment.  If  the  temperature  exponent  is  set  equal  to  zero,  it  is  clear 
that  a  much  stronger  feature  is  predicted.  Choosing  a  temperature  exponent 
between  zero  and  0.5  would  probably  resolve  the  discrepancy,  but  this  would 
be  physically  unrealistic.  This  temperature  dependence  would  then  become 
Inconsistent  with  what  P-branch  data  are  available,  and  does  not  seem  theo¬ 
retically  justifiable.  The  j-dependence  of  the  llnewidths  is  so  weak  at 
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high  temperature  that  this  does  not  seem  to  be  a  sensitive  parameter.  As 
will  be  seen,  collisional  narrowing  calculations  at  higher  pressure  are  in 
reasonable  agreement  with  experiment,  and  thus  it  does  not  seem  likely  that 
the  fault  lies  with  the  collisional  narrowing  models.  For  a  given  rate  law 
model,  the  question  of  the  proper  temperature  dependence  loses  importance 
as  the  pressure  is  increased  because  the  calculations  lose  sensitivity  to 
isolated  linewidth. 

I  One  possibility  that  should  be  considered  is  the  linewidth  augmentation 

effect  suggested  by  St.  Peters  (Ref. 32).  He  argues  that  the  Raman  line- 
widths  are  effectively  increased  by  the  sum  of  the  finite  axial  mode  widths 
in  the  pump  and  Stokes  lasers.  That  is,  he  asserts  that  the  effective 
linewidth  which  goes  into  the  third  order  susceptibility  calculation  is 
r  =  I'j  +  Tp  +  Fg,  assuming  that  all  of  the  modes  in  the  broadband  pump 
I  and  Stokes  have  roughly  the  same  widths.  This  argument  is  tempting  because 

there  is  speculation  that  the  Stokes  modal  widths  might  be  as  much  as 
.05  cm~^,  which  is  the  kind  of  number  that  is  needed  to  resolve  the  dis¬ 
crepancies  discussed  above.  This  explanation  is  also  consistent  with  the 
data  at  high  pressure,  and  if  it  were  true  would  have  very  significant  im¬ 
plications  because  the  linewidth  at  low  pressure  would  then  become  an  in¬ 
strumental  feature  and  not  a  function  of  medium  conditions.  This  argument 
is  being  raised  here  only  as  a  possibility.  It  is  questionable  however, 
because  it  does  not  need  to  be  invoked  in  the  case  of  N2  spectra  obtained 
from  the  same  apparatus  where  the  observed  interference  at  the  bandhead  is 
,  significantly  smaller. 

I 

Another  possibility  concerns  the  longitudinal  temperature  gradients 
which  are  known  to  exist  at  high  temperatures  (Fig.  6).  In  propagation 
through  such  gradients,  each  Fourier  component  of  the  anti-Stokes  field 
will  be  given  by  a  coherent  superposition  of  the  contributions  from  a 
I  volume  at  each  temperature  along  the  path;  thus,  interference  effects  will 

occur  that  might  be  exaggerated  by  the  different  temperature  weightings 
assigned  to  the  susceptibily  at  each  location.  The  bandhead  has  been  seen 
to  be  quite  sensitive  to  such  interference  effects.  However,  the  problem 
with  this  argument  is  that  the  measured  gradients  are  not  severe,  as  can  be 
seen  in  Fig.  6.  Certainly  no  ambient  CO  contributions  are  being  mixed  in, 
nor  is  the  bandhead  spike  predicted  at  the  lowest  temperature  prevailing  in 
the  interaction  region,  and  an  incoherent  superposition  of  spectra  with 
temperatures  within  the  interaction  length  do  not  show  the  interference  at 
the  bandhead.  There  is  the  problem  too  of  reconciling  this  with  the  good 
fits  obtained  at  higher  pressure.  Although  the  gradients  might  become  smaller 
temperature  distributions  at  high  pressure  have  yet  to  be  determined. 

Thus,  the  interference  feature  which  is  observed  at  one  atmosphere  in 
the  hot  spectra  is  still  something  that  is  not  entirely  understood,  and 
there  does  not  seem  to  be  a  completely  satisfactory  explanation  for  it  at 
the  present  time.  In  any  event,  the  failure  of  the  theory  to  account  for 
it  does  not  necessarily  have  negative  implications  for  diagnostics.  The 
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theory  and  experiment  agree  well  in  all  other  parts  of  the  spectrum,  so 
that  regression  analyses  will  yield  accurate  temperature  predictions,  and 
measurements  that  rely  on  the  total  signal  strength  (l.e.  the  integrated 
area  of  the  spectrum),  will  be  little  affected  by  the  narrow  interference 
feature  at  the  bandhead. 

The  pressure  dependence  of  the  295K  spectra  has  previously  been  com¬ 
pared  to  the  predictions  of  the  SPL  and  the  Gordon  model  (Ref. 17).  These 
results  are  shown  in  Fig.  10.  Both  these  models  require  an  average  elastic 
dephaslng  contribution  of  about  5%  of  each  isolated  linewldth  to  be  consis¬ 
tent  with  the  data  at  high  pressures.  In  reality,  it  is  likely  that  the 
elastic  contributions  would  have  a  j-dependence ,  like  the  isolated  line- 
widths.  Based  on  the  results  of  Koszykowski  (Ref.  33)  it  would  be  reason¬ 
able  to  expect  that  the  elastic  contributions  would  be  most  important  for 
the  low  j-values,  with  diminishing  importance  as  j  increases.  The  calcu¬ 
lated  results  for  the  EPL  model  at  one  hundred  atmospheres  are  shown  in 
Fig.  11.  As  can  be  seen,  the  agreement  is  impressive,  although  it  too 
would  require  some  pure  broadening  (elastic)  contributions  to  make  the  fit 
perfect  (the  experimental  bandwidth  is  1.58  cm”^ ,  the  theoretical  band- 
wldtli  is  1.35  cm~^). 

In  Figs.  12-20  theory-experiment  comparisons  are  presented  for  condi¬ 
tions  which  represent  the  range  of  the  experimental  investigations.  The¬ 
oretical  overlays  corresponding  to  both  the  SPL  and  EPL  are  presented;  for 
the  hot  calculations  at  modest  pressures,  the  perturbation  expansion  Eq.  6 
was  employed;  for  pressures  in  the  vicinity  of  100  atmospheres,  the  full 
G-matrix  inversion  was  necessary.  There  may  be  a  general  tendency  for  the 
EPL  model  to  overpredict  the  extent  of  narrowing;  the  exception  is  Fig.  20, 
where  it  seems  to  give  a  better  fit  than  the  SPL.  On  the  other  hand,  there 
are  certain  spectra  where  the  SPL  does  not  predict  quite  enough  narrowing. 
Nevertheless,  these  defects  are  small,  particularly  when  one  considers  the 
enormous  errors  that  would  be  incurred  by  neglect  of  narrowing  entirely 
(Figs.  18,20).  It  seems  reasonable  to  conclude  that  both  models  do  quite 
well.  The  Gordon  model  predictions  are,  in  general,  slightly  less  satis¬ 
factory  than  those  of  the  simple  power  law.  As  an  example.  Fig. 16  compares 
a  prediction  of  the  Gordon  model  with  those  of  the  power  laws  at  lOA  atmo¬ 
spheres  and  645  K.  Elastic  contributions  were  not  included  in  any  of  the 
hot  spectra;  the  hot  N2  calculations  discussed  earlier  indicated  no  need 
for  such  contributions  in  that  molecule.  As  discussed  earlier,  the  ques¬ 
tion  of  the  proper  temperature  dependence  of  the  linewidths  becomes  moot  at 
higher  pressures;  numerical  tests  for  the  hot,  high  pressure  spectra  (Figs. 
18,  20)  showed  no  significant  sensitivity  to  linewldth  for  each  model. 

The  results  of  these  comparisons  suggest  that  it  might  be  desirable  to 
consider  a  rate  law  which  combines  the  SPL  and  EPL  in  some  way.  The  prob¬ 
lems  in  the  fundamental  bandhead  of  hot,  low  pressure  spectra  also  are  a 
question  that  should  be  addressed.  It  might  be  useful  in  this  regard  to 
undertake  ab  initio  calculations  of  the  kind  that  have  been  performed  for 
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N2  (Ref.  27).  The  investigations  undertaken  in  this  contract  represent 
quite  an  extensive  test  of  the  basic  theory  of  high  pressure  CARS,  and  it 
is  encouraging  that  such  good  results  can  be  obtained  from  llnewidth  models 
of  such  simple  form.  The  diagnostic  Implications  of  the  work  completed 
under  the  current  contract  are  encouraging.  Based  on  the  results  of  these 
pure  CO  investigations,  it  is  expected  that  minority  species  CO  spectra 
will  also  be  successfully  Interpreted  by  the  collisional  narrowing  model. 
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AREAS  OF  FUTURE  WORK 

The  successful  completion  of  the  tasks  In  this  current  contract  have 
resulted  In  verification  of  the  ability  of  the  UTRC  CARS  computer  code  to 
accurately  model  and  predict  the  CARS  spectrum  of  carbon  monoxide  over  the 
temperature  range  300K-1500K  and  for  pressures  from  1-100  atm.  The  code 
now  represents  a  powerful  tool  for  accurate  temperature  and  concentration 
measurements  of  this  molecule  In  combustion  media.  A  logical  extentlon  of 
the  current  work  Is  to  proceed  with  experiments  In  high  pressure  laboratory 
flames,  (l.e.  CO/alr)  and  attempt  flame  chemistry  studies  through  CARS 
measurements  of  temperature  and  species  concentration  In  high  pressure 
flame  environments. 
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EXPERIMENTAL  CARBON  MONOXIDE  CARS  SPECTRA 
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paBHaBBBI 

liaOBBBBBI _  .  -  _ 

kiaHBBBBBBBBBBBBBBBaaaBaaaBBBBBBBBBI 
paBBBBBBBBBBBBBBBBBaaaBaaaBBBBBBBBI 
paUBBBBBBBBBBBBBBBBaaiiaaBBaBBBBBBBI 
paaaaBBBBBBBBBBBaBaaaapaaaaBBBBBBBi 
NaaBBBBBBBBBBBBaaiiHasisaaaBaBBi 
- IHBI 


IBBk.^= 


"31 


..laaBflBBBBBBBBBflBBBaaaaaaBaBHBBilBB _  _ 

paaflaBBBBBBBBBBBBHaaBBaBaBBBBaiSBBBB^rr^^^BiflBl 
laBBBBBBBBBBBBBBBBBBHaaBaaBBBBBBBBBl^^  HiiBBBgii 
laaaBBBBBBBBBBBflBHBaaaBaaBffiiiBBBBBE^^^BBnBBn 
bj^aaaaaBBBBBBBaBBaaBs^aaeaaBBBaBBr^aBBBiiiiBI 

paaBaBBBBBBBflBBBBBaBaaaBRIBBBBaBfe.  _ 

paaaaasaaflBBBflBBeaaai!PPiiaaBBBaBiKraBBBBBBBaBi 
l^iiiiiiiBaBBBaBBBaBaBBaBiiiiiiaHaBaBBBr'-^aaBBBBBir' 
hH::llia«BBBBBBBBBBBaBBiiaiSBBBBIIB»ii!!’:::;:.«BBaBBBBBil 
paaiiaaBBBBBBflBBBB'-’-T-^^.i.iMaBBBBBBflaBBaBaBBI 

paaBBBBBBBBaHBBiiaBaaapaaBBaBBBaBBBBBBaBBBaite^llij 
paBBHBBflBBBBBBBHBaaHaaaaBBliBBBBBBiBBBBBBBBeBB^Bi 
hiaaBBBBBBBBBBBBBBBaBHKaaHBBBBBliflBBBBBBBaBaBini 
paaBBBBIiaBaBBBBBBBiga^aaKBliiailllBaBBBBBnBBBBBiiBkl 
ps^^^^BBBBBBBBBBBliBBaiiBaBBIiaBBBBBBBr 
Jliaiii|IBflBBBBBBBBBBHflBaHBBBaBitlHBBBBB~~ 
|u^aaBBBBBBBBBBBBBBBBBBBBBBBnBBBBBB„, 
paaBBBBBBBBBBBBBBBBBHBBSiaiiinaBBBBBBI 
paBBaBBBBBBBBBBBBBaaBHaijiaaHBaBBBBBBBBflBaBBBilMBII 
pBBaBBBBBBBBBBBBaBBSieiaBaBHBBnaBBBBBBBBBBBBBiBBBii 
laeaBBBBBBBBBBBBBBaaiiiiBBaaBnBBBBBiMIBBBBI 
pai^aaBBBBBBBBBBBiilSa^aaiiKHiaaBKaaBBBBBBBBI 
Igi^jjS^aBMISBBBBBBaBBHBBaaaMIBBBBBBBBBBiMiBBI 
NHBaaaBaBBBBBBBaaBaBiiBBBMaBBBBaBBBBBpwi 
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laaapHaaBKaisaaaiiHBHaaHBiiiiaHaBaHiigaaaasBaBaHHaHaSq 

^aigiliiiBi^aaaaaaaaaaisaaaaaaaaaaaaaaaBaaw 
^^^■■■■■■naaBBaaaaa^iiaaaiiiaaiiaaaaiiaaBar^- 
,  iBaaaaaBigaiiissaiiiBsgaaaaaaaayoiMags 

la^^^^^B^KaBaaaaaBaaaaaaaaaaaaaaaaaaaBeaaaff’ 

aaaaBaBBaaBaHaBBBaaaaaaaaBacawji 

,  HMBaHHBB— Bwaamaia 

fWlTT""""""”""*""^*"*""*"*""”""**" 

la^^^^^H^Haaaaaaaaaaiiaaaaaai - - - 

^M^^^EmaaaaaaaaaaaaaaaaBi  _ 

^ _ aaaaBBBBiaaiaSSBBBBBSBBBaBEa! 

issiissassssssssEssssgsassB 

aflBSiBBBflBBBBBBBHBBBBBBBBaBBaiHBBBBBBIlIBBBBBBBBB'tfl 
SBBBBBBBBBBBBBBBBBBBBBgliliaBiiliBBBBBBBBBBBBBBBBBi:: 
aaaHilBBaBBBBBBBBiaHBflBBBaBBilBBBBBBBBHBBBBBBBBBBI 
IBBBBBBBBBBBBBBBBBflBBilBaBBBBBBBBBBBBBBBBI 
aagBBBBBaBBBBBBBBBBBaa^BaBaBBBBBBHIiliBBHBI 
liiaBaBSiBBBBiiBBflflBBflBflBBiiBiiBiillilBBIlBBBBBflBBBL 
sSaBaBBBBBigaBBggBBBaBaaiBBafelliBiaBiiBBBBHBBBBBBBBBBM 
^BaBaaBaaaBeaaBBBaBBaiiaaBjiiBaiaaiaaBaaaBaaaaBiH^aii 
IBflBBBaBflBflBflBBBBBBBiliiBiiaBllillBBBBBBBBBBBBBBBiiBr^r^^ 
aBBliBBBBBBBBBBBBBBBBBBBliBBBBBBflflBBBBBiBSBanBL'll 
liilBBBBBBBBBBBBBBBflBBBBBBBBBBBBBBMBBBBBfll 
^BBBBBBBBBBBBBBBBBBBBBIBBaBBBBBBBBBBBBBBfll 

BUBgBBBflBBBBBflBBBflBBflBBBBBBBBBBBBaflBgil - 

SSBBBliBBBBBflBBBflBBBBBBBBiieBBBBBflBBBBMI 

aBBBBliBBBBBBBBflBBBflBHBBb^BBBBBBBBI - 

BBBBBBBBBBBBBBBBBBBIiliaiBBBBBBBBBBI 
gaBBBBBBBBBBBBBBBBflBBBBHaiiBBBaBBI 
aBUBBBBBBBBBBBBBBBBBBflBflBBBBBflBBi 
sBBBBBBBBBBBBBBflBBBBHBBBBBBBBBBBI 

aayjBBBBBflBBBBBBflBflBBBBBBBBBBBflBfll _ 

HaBBliBilBBBBBBBBBflBaBaBBliiBBBBBBBBBBBBBBBBBa^UMi 
JB^BBaBBBBBBBBBBBflBBBBBBBBBilBBBBBBflBBBBBa»:^BiBl 
ii  iniininim  TPiimiiMiiiBBBfiiBTW  iii  iii  iii  r  numip i  hiiiiiibbbbiibbbbbiibi  bbbBi 

aBBBBBBBBBBBBBaiBBBBiiBliiiBiiBii3i!!BaBBBBflBBBBBB:3<«BBBBl 

^^nBBBBBBBBBBflBBBBBBBBBiiBiiBBiiaiBBBBBflBflBBBa^ - 

IBBBBBBBBBBBBBBBflBBfliiBHBBBBflBBBBBBBBBBBB’^ _ 

H|llgBBBBflBBBBBBflBflBiaBBBBBBBaBBBBBBBaiBBer^BBI 
iB!iBBBBBBBBflBBBflBBSBBBaiBBBIIiBHBBaBBBiii&r;.«Bfl9 

l:lll3iB»BBBBBBBBBflBBBBBBBBBiiBBBBBBBBBBI - 

-JiiSiiiiaiSBBBBBBBaBBBBBBBBBaBBBBBBBBBBBBI _ 

iniilBBBgBBBaBBBaBBBBB^isPJBIlBBBBBBBBBBBr^JBI 

liB9BBiiiiifBBflBBBBBBBBaBBBBBBBBBBBBBBP»:L«aBBI _ 

iaBpiBBBBBaBBaaBBBBBBaBBaBBBaaagr^MMiaaiWBaaBBBi 
32!rJifBaiiBBBBBBBBaBBBaaaBBailBBaP'*”^ciaBBHlBBBBBaaBBfll 

i::aiiBaaBBBaaBaBBBiasaHaaBBBBssi°":^».aBaBBBi - 

iiBBB9aiiaBBBP«*"!;=:===;^;;^asaiiBnBBBBBBBBBBI 
SailBBBBBBBBBIi^Bn»«nB»n«^=Z:=:=;;P»!««Bai 
I^BaBaBBBBBBBBaBBBBBgBaaBBBBiiBBI 
,  iBBaaaBBBBBBBBBBBBlIBBBaiaaiBiiiai— 
l^i:lBil.aBBBBBBBBBBBBBaa^^aBaaBBilBaBBBBBBaBBBraBiik1i 

,  ii^BaBBBBBBBBBBBaBBHBiiBBBaBBBa - 

L:^^ae£cuauHMH|flHBBBBaMBBBBBBaBI 

.aa,«»».jaaa8SggiiisssgBBaBBBB88i« 

iiaaaaaaaaEaiaaaaaaaaaaaassai 
iiisaaaaaaaaisaiiiigagaaaaaaaair 

liiaaaaaaaasiaaap"™— « 


A-17 


lasBBHBSsiasiunni 


BBBBaHBBBSBHaBBBBaw^ 


_ IBBBBBHBHIiBilBBBBBI _ 

IBBBniBBBBBBBBBBliBBBBBBBBBBBBB: 

_  _ JBBBBBBBBBBBBBBHBBBBBBBBBBBBB' _ 

IBBBBBaBSaSBSSSBaaBSBBBBBE; 

■BBBBBBBBBBBBBBBIHBaBBBiBBBBBBI 

IBBBBiBBBBBBBflBHBBBflBBBflBBBBl _  _ 

IBBBBBBBBBBBBBBHBBBBBBHBBBBBBnillM^^I 
■"■BUBBnBBBBHBBBBflBBflnBflBiSE-^BF-' 
IBBBBBBBBBBBBBBBBBBBBBBBBPWT-r^l 
IBBiBiBBflBBBB»BBBBBBBBBBUiBB»2=^l 


IBBBBBBI 
IBBBBBBHI 

Ibbbbbbbbi 

IflBBBBBBBI 


■■aBBI 


IBBBBI 

Ibbbbi 


IBBBBI 

Ibbbbi 


IBBBBBBBBBBBBBBBBBBBBliBrEE^^-MMBBBBBBinBBI _ , 

'BBBBBBBBBBBBBBBBBaiBBV»^===ii«BBBnBBBBBBMBB 

BBBBBBBBBBBBBflBBBB»**»==;;;i.«aBniBBBBBBBBBBB 

BBBBBBBBflBflBiiB!>Ew««!!--^^»«BBBBBBBBBB»nBB 


IBBI 


;izi:rr:E=~«iaiBBBBBi 


EBBBBBB^  .IBBBBBBI 

■^^fS^aBBflBBBBBBI 


iBBaa 
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lawBiinHi 
lisaBBiBr 
IBBBBBBI 


lanRSBBB 
IHBSSSSREBSIHSSSaaS 

IBBBBBBBBBBBHBBBBBBBBBBBBBBBBBBBtr 
IBBBBBBBBBBBBBBBBBBBBBiiBiBiiBBBB 

^  IBBBB  “  "  ^ 

IBBBB 

1B8R 

IBBBfllBBBBBBBBBBBBBUBBBflBBI 
IBBBBBBBBBBBBBBBBBBiiBiBiiBI 
IBBBBBBBBBBRBBRtBBBBBBBBBBBBI 
IBBBBBBBBBBBBBBBBBBBBBBBBBBL-. 

IBBBBBBBBBBBBBBHMBBBBBBBBBBBBi _ 

IBBBiBBBBBBBBBgBBSBBBBBiiBBiBeBr«i 
IBBBBjjBBBBBBBHHBiiBBBBBBBBBBHBlfc'Ril 

_ IBBBBBBBBBBBBBBBBBBBBBBBBBBflr^flBl 

IBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBB^l 
IBBBBBBBBBBBBHBBBBSBBBBBBBBBBBBHBBB.' 
IBBBBBBBBBBBBBBBBBBBBBBBBBBr 


laBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBlBBBBBBBBk:' 

lBBBBBBBBBBBBBBBBBBBBiBBBBBBBBBBflBBBBBBBBBBBr^_ 

|ffiSBBBBflEBBBBBBBBBBBBBBBBBBBBBiiBBBBH!lBBBBBBBP*^^l 

|ailBBBBBBaBBRBBBBBBBBBSBBBBBaHBBaBBBBHBflBi9BBb.;rBI 

BaSBSSSBSBSRBSSIBiEBBBBBHRig'  ~ 

BBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBB^' 
BBBBflBBBBBBBBBRBBBBBBBBBBBBBBiiBBB;^^« 
BBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBB 

_ BBBHBBBBBBBBBBBBBBBBBBBBBBBBBBBBI 

BBBBBBBBBBBBBBBBB9BBBBBBBBBBSBBBBBBHBBBBI.  ^ 
giaBBBBBB9BBBBBBBBBBBBBBBBBBRBBBBBBBBBBBB«;r«Bi 
I^SBBSIBSBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBHCr-irtil 
iBBBBBBBBBBBBBBBBBBBBBBBHHBBBBBBBBBBBBBBK'^^dBBBI 
BBaBBiBBBBBBBaBBBBBBBBBBBBBBBBBBB*~~ 

I— — BBBBBMBBWBBBiBBBBBBBBBBBiaBBBii»> 
iBBBBBBBBBBBBBBBBBBiBBBBBBBBBBBBBBBBBBifcg 
iBBBBBBBBBBBBIBBBBBBBBBBBBBBBBBBBBBBBBBi^, 
iBBBBBBBBBBBBBBBBBBBBBBBBBBBBBflBBBBBflBk-'T! 
iBBBBBRBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBkn 
iBBBBBBBBBBBBBBBBBBBBBBBHBBBBBBBBiBBBF'^BBBiliB 
iBBBBBBBBBBBBBflBBBBBBBBBRBBBBBBBBiBBBr^  MBMBI 
iBBBBBBBBBBBBBBiBBBHiBBBBBBBBBBBBBBir.^! 
iBBBBBBBBBBBBBBiiniiiBBBBBBBBBBBBBP' 
iBBBBBBBBBBBBBBiiBiBBiBBBRBBBBBBBS^'-^BBj 

l8PHBBBBBB«iBBBBBiiBBBiBBBBiBBBir~:;Sii.aiBBBi _ 

|gBBBBBBBBBBBBBBBBBBBBBBBI»":i^ik«BBBBBBBBBBBBBBI 
lBBBBBBBBBBBBBBBBnB*«!^!l..aBBBflBBBfBBBBBBBBBBI 


iBBBaaBaai 


iBaSBI 

IRBBai 


laaaai 

Ibbbbi 


iBBaBBaaBaaaaaaB 
iBBBaaaiiiiaiBBB 
BaaaaaBBBi 

aaflaaaBBBi 

BEEBBSBBI 
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